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{Sj ■ Abstract 

■ It is shown that every d-by-d unitary matrix can be written as the product of two-level unitary matrices 



with additional structure and prescribed determinants. The result is then applied to show that every quantum 
(unitary) gate acting on n-qubit register can be expressed as no more than 2 n (2™ — l)/2 controlled single 
qubit gates chosen from 2™ — 1 classes so that the quantum gates in each class share the same n — 1 control 
qubits. Related results are discussed. 
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1 Introduction 

The foundation of quantum computation [9] involves the encoding of computational tasks into the 
temporal evolution of a quantum system. A register of qubits, identical two-state quantum systems, 
is employed, and quantum algorithms can be described by unitary transformations and projective 



d 



measurements acting on the 2 n -dimensional state vector of the register. In this context, unitary 



transformations (matrices) are also called quantum gates. 



It is well known [1] that the set of single qubit gates and CNOT gate are universal. In other 
words, any unitary gate acting on an n-qubit register can be implemented with single qubit gates 
£NJ ! and CNOT gates. Typically, the proof can be done in the following steps; for example, see [8, §4.6] 
and [9, pp. 188-193]. 

Step 1. Every d x d (special) unitary matrix can be written as a product of no more than 
^ ' d(d — l)/2 two-level (special) unitary matrices. (A two-level unitary matrix is obtained from by 
changing a two-by-two principal submatrix.) 

Step 2. Every quantum gate of size 2 n corresponding to a two-level unitary matrix U can 
be converted to a controlled qubit gate U using at most n—1 swaps of states by the Gray code 
technique. (Here each swap is a product of single qubit gates and CNOT gates.) Then the controlled 
qubit gate U can be expressed as the product of single qubit gates and CNOT gates. One can then 
reverse the swaps of states to get the decomposition of the unitary U into single qubit gates and 
CNOT gates. 

There has been considerable interest in finding efficient and practical algorithms in constructing 
unitary gates; see [2, 3, 5, 6, 7, 8, 9, 12] and the references therein. The purpose of this note is to 
show that one can express any unitary gate (matrix) acting on n-qubit register as no more than 
2 n (2 n — l)/2 controlled single qubit gates. These gates can be chosen from 2 n — 1 classes such that 
each of these classes consists of controlled qubit gates using the same n — 1 control qubits. 



1 Corresponding author. 
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While our proposed construction does not improve much in terms of the complexity of the 
decomposition of unitary gates (because controlled qubit gates are still expensive to implement), 
it gives a more elegant proof for the universality of (selected) controlled qubit gates. Moreover, 
our proof uses only elementary properties of unitary matrices, and a clever combinatorial idea to 
decompose a given unitary gate as the product of (selected) controlled qubit gates. The method 
is easy to understand and implement. For example, a Matlab program for decomposing a 4-by-4 
unitary matrix into no more than 6 controlled qubit gates by Ms. Rebecca Roberts is available at 
http: / /cklixx. people. wm.edu/mathlib. html. 

To make our method more portable to other problems, we present a general lemma which 
allows one to decompose a unitary matrix to two-level unitary matrices with additional structure 
and prescribed determinants. The technique and insight are useful in other problems in quantum 
information science and related subjects. 

We will present our results, proofs, and examples in Section 2, and give a short discussion in 
Section 3. 

2 Results, Examples, and Proofs 

Let P = {31,32, ■ ■ ■ ,3d) be such that the entries of P correspond to a permutation of (1, 2, . . . , d). 
A two- level unitary matrix is called a P- unitary matrix of type k for k £ {1,2,... , d — 1} if it is 
obtained from 1^ by changing a principal submatrix with row and column indexes jk an d jk+i- F° r 
example, if P = {ji, 32,33,34) = (1, 2, 4, 3), then the three types of P-unitary matrices of type 1, 2, 
and 3 have the forms 

/* * 0\ 
* * 
10 

Vo 1/ 

respectively . We have the following. 

Lemma 2.1 Let P = {31,32, ■ ■ ■ ,3d) be such that the entries of P correspond to a permutation of 
(1,2,... , d) . Then every d-by-d unitary matrix U can be written as a product of no more than 
d(d — l)/2 P -unitary matrices. Moreover, those P -unitary matrices can be chosen to have any 
determinants with modulus 1 as long as their product equals det(U). 

An immediate consequence of the lemma is the following. 

Corollary 2.2 Every d x d special unitary matrix can be written as a product of no more than 
d(d — l)/2 P -unitary matrices each of which has determinant 1. 

It is instructive to illustrate a special case of the lemma. We consider the case when d = 4 
and P = {31,32,33,34) = (1,2,4,3) as above. (As we will see, this example is relevant to our later 
discussion on quantum gates.) 
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be a four-by-our unitary matrix. Let //i, /X2, fJ,e be such that fii, fj,2, p-e € : |z| = 1} and 
■ ■ ■ A*6 = det(CJ). We divide the construction into two steps. 

Step 1. We consider the column of U labeled by the first entry of P (i.e., the first column). 

Choose P-unitary matrix U\ of type 3 with det(£7i) = fii such that the (j'4, ji) = (3, 1) entry of 
U\U is as follows. Let «i = \/l a 3i | 2 + [041 1 2 and 
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Next choose P-unitary matrix U2 of type 2 such that (j'3, ji) 
follows. Let u 2 = 



u 32 u 33 u 34 
a 42 a 43 a 44 ' 



(4,1) entry of U 2 U 1 U is as 



0-21 1 2 + u\i and 
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Then U 2 U 1 U = 
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Now, choose P-unitary matrix U3 of type 1 so that the (j'2, ji) 
also 0, and the (1, 1) entry of U3U2U1U equals 1 as follows. Let 



(2,1) entry of U 3 U 2 UiU is 
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Then V = U3U2U1U 
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Note that the first row of V has the form (1, 0, 0, 0) because V is unitary. 

Step 2. We turn to columns of V labeled by j'2 = 2 and j'3 = 4. 

Choose P-unitary matrices U4, C/5 of types 3 and 2 with determinants p, 4 and ^5, respectively, 
so that the (34,32) = (3,2) entry of U4V and the (33,32) = (4,2) entry of U5U4V are 0. Then 
choose a P-unitary matrix U§ of type 3 with determinant p,$ so that the (34,33) = (3,4) entry of 
UQU5U4V is 0. Here are the zero patterns of the matrices in the process: 



U4V = 




, u 5 u 4 v = 



1 0^ 

10 

0* 

,0 0* 



, U 6 U 5 U 4 V = 



1 0* 

10 

10 

,0 1, 



Note that in the last step, the (j'4,.7'4) = (3, 3) entry of UQU5U4V is 1 because 
det([/ 6 [/ 5 W3IW1) det(£/) = £ 6 /2 5 • • • £1 det(l7) = | det(?7)| 2 

Consequently, 



1. 



Clearly, each l/j is a P-unitary matrix of the same type as Uj,j = 1, 2, 6. 

Obviously, we can skip some of the P-unitary matrices if the entry to be eliminated is already 
during the process. We can now present the proof of Lemma 2.1. 
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Proof of Lemma 2.1. Let P = (31,32, ■■■ ,3d) where the entries of P is a permutation of 
(1,2,... , d) . Let U be a d-by-d unitary matrix. We extend the construction in the example to 
the general case as follows. 

Step 1. First consider the j\th column of U. One can choose P-unitary matrices Ud-i, Ud-2, ■ ■ -, 
XJ\ of types d — l,d — 2, . . . , 1 with prescribed determinants to eliminate the entries of U in po- 
sitions (jd,ji), (jd-i,3i), ■ ■ ■ , (j2,ji) successively, so that the jith column of the matrix V = 
U1U2 ■ ■ ■ Ud-iU equals the jith column of Id- Since V is unitary, the jith row of V will equal 
the jith row of Id- 
Step 2. Consider the ,72th column of V. One can choose P-unitary matrices V^-i, Vd-2, . . . , V2 
of types d— l,d — 2, ...,2 with prescribed determinants to eliminate the entries of V in positions 
(jd,h), {jd-1,32), ■ (33,32) successively, so that the j 2 th column of the matrix W = V 2 ■ ■ ■ V d -iV 
equals the ,72th column of Id- Since W is unitary, the ,72th row of W will equal the ,72th row of Id- 

We can repeat this process in converting the ,73th . . . , jd-ith columns to the ,73th . . . , jd-ith 
column of Id successively, by P-unitary matrices. Finally, the (jd,3d) element will also be 1 by 
the determinant condition imposed on the P-unitary matrices. Multiplying the inverses of the 
P-unitary matrices in the appropriate orders, we see that U is a product of P-unitary matrices as 

asserted, and the number of P-unitary matrices used is no more than (d — 1) H h 1 = d(d — l)/2 

because some of the P-unitary matrices may be chosen to be identity if the entry to be eliminated 
is already during the process. □ 

We will apply Lemma 2.1 to the decomposition of a quantum gate into the product of controlled 
qubit gates. In this application, we need a special selection of the vector P = (i\, 12, - - - , id)- 

In the following discussion, we always assume that n is a positive integer and N = 2 n . Recall 
that a Gray code G n [10] is defined as an iV-tuple G n = (X±, X2, ■ ■ ■ , Xn) such that 

(a) each X\, X2, ■ ■ ■ , Xn are length n binary sequences corresponding to binary representation of 
the numbers 0, 1, . . . , N — 1, arranged in a certain order, 

(b) two adjacent sequences Xj and X/+i differ in only one position for each j = 1, 2, . . . , N — 1, 
and 

(c) the sequences Xn and X\ differ in only one position. 

One can construct a Gray code G n recursively as follows. 
Set Gi = (0,1); for n > 1 and N = T, if G n = (X lt . . . , X N - U X N ), set 

G n+ \ = (OXi, . . . , 0Xjv-i, OXn, IXn, 1Xjv-i, . . . , lXi). 

For example, we have 

G 2 = (00, 01, 11, 10), G 3 = (000, 001, 011, 010, 110, 111, 101, 100), etc. 

Let us adapt the definition of P-unitary matrices to define G n -unitary matrices, which corre- 
spond to controlled single qubit gates in quantum information science. To this end, label the rows 
and columns of an N x N matrix by the binary numbers 

0---0, 0...01,..., 1---1. 
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An N x N two-level unitary matrix is a G ra -unitary matrix of type k if it differs from In by a 
principal submatrix with rows and columns labeled by two consecutive terms X\~ and X^+i in the 
Gray code G n = (X\, X2, Xn), k £ {1, 2, . . . , N— 1}. Clearly, there are N—l types of G ra -unitary 
matrices. Since X^ and X^+i differ in only one position, every G n -unitary matrix corresponds to 
a controlled single qubit gate. 

For example, (SVunitary matrices of types 1,2,3 have the following forms, respectively: 

00 01 10 11 00 01 10 11 00 01 10 11 
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Note that these are exactly the P-unitary matrices when P = (1, 2, 4, 3). 
It is now easy to adapt Lemma 2.1 to prove the following. 

Theorem 2.3 Let n be a positive integer and N = 2 n . Every N x N unitary matrix U is a product 
of m G n -unitary matrices with m < N(N — l)/2. Furthermore, the G n -unitary matrices can be 
chosen to have any determinant with modulus 1 as long as their product equals det(£7). 

Proof. Identify (0, 1, ... , N—l) with the iV-tuple of binary numbers (0 • • • 0, • • • 01, . . . , 1 • • • 1); 
label the rows and columns of U by the binary numbers • • • 0, • • • 01, . . . , 1 • • • 1. Then apply 
Lemma 2.1 to U with P replaced by G n . □ 

Since every G n -unitary matrix corresponds to a controlled single qubit gate, by Theorem 2.3 
every unitary gate can be written as product of controlled single qubit gates from 2 n — 1 classes so 
that matrices in each class share the same n-1 control qubits. 

The example before the proof of Lemma 2.1 can be viewed as an illustration of Theorem 2.3 when 
n = 2 if we identify P = (1, 2, 4, 3) and G2 = (00, 01, 11, 10). An illustration of our construction for 
n = 3 is given in the Appendix. 

3 Discussion 

In this note, we obtained a decomposition of a d-by-d unitary matrix as product of special structures 
specified by a vector P = j'2, . . . ,jd) such that the entries of P correspond to a permutation 
of (1,2, ... ,d). The result is then applied to show that every unitary gate U acting on n qubits 
can be decomposed as product of special two-level unitary matrices corresponding to controlled 
single qubit gates. This was done by using Gray code G n = (X±, X2, ■ ■ ■ ,Xn) with iV = 2 n , and 
constructing G n -unitary matrices which are two-level matrices obtained from In by changing its 
principal submatrix with row and column indexes X^ and Xf, + \ for k = 1,2, . . . ,N — 1. 

There are other applications of Lemma 2.1. For example, if P = (1, 2, . . . , d), then P-unitary 
matrices are two-level tridiagonal unitary matrices. In numerical linear algebra and other appli- 
cations, it is useful to decompose a matrix into tridiagonal forms with simple structure; e.g., see 
[4, 11] and their references. In our study, the focus is on quantum information science, and the 
decomposition of the unitary matrix respects the tensor (Kronecker) product structure of matrices. 
The same technique may be useful for decomposition of matrices with other multilinear structures. 
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Appendix: Illustration of Theorem 2.3 for three qubit unitary gates 



Let us illustrate Theorem 2.3 when n = 3 in the following. It will further demonstrate how to 
eliminate the nonzero entries of the given unitary matrix using the G n vector, which will deter- 
mine the order of the columns to be treated, and the order of the entries to be eliminated in a 
selected column. Evidently, depending on the applications, one may use different Gray codes for 
the decomposition of a given unitary gate or a special class of unitary gates. 

To illustrate Theorem 2.3 when n = 3, note that 

G 3 = (X U X 2 , X 8 ) = (000, 001, 011, 010, 110, 111, 101, 100). 

Suppose U is an eight-by-eight unitary matrix. We divide the decomposition of U in several steps. 

Step 1. First, we consider column X\ = 000 (i.e., the first column). Choose G^-unitary matrices 
U^l, Uxl-i Uxl-i Uxl, Uxl-i Uxl-i where the upper index in each symbol corresponds to 

the type of the GVunitary matrix, with appropriate determinants (as done in the case when n = 
2) to eliminate the nonzero entries in column X\ with row indexes Xg, Xj, Xq, X§, X4, X3, X2, 
successively. Here note that if we identify the row and column indexes (000, 001, . . . , 111) with 
(1,2,..., 8), then the Gray code (X±, X2, X3, X4, X5, Xq, X-j, X$) corresponds to (1, 2, 4, 3, 7, 8, 6, 5). 
So, we are eliminating the entries in the first column lying in rows 5,6,8,7,3,4,2, successively. Let 
U Xl = UgU%UgU%UgugU* Then 



u Xl u = 



000 001 
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Step 2. We consider column X2 = 001 (i.e., the second column). Choose GVunitary matrices 
Uxl , Uxl , Uxl , Ux^ > Uxl , Uxl with appropriate determinants to eliminate the nonzero entries in 
column X2 with row indexes Xg, X7, Xq, X§, X4, X3, successively. (Here, using the usual row and 
column indexes of an eight-by-eight matrix, we have eliminated the nonzero entries in the second 
column with row indexes 5,6,8,7,3,4 successively.) Let U Xi = ^x^xt U xt U ^ U xl U xl- Then 



U X2 U Xl U 



(h 
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* * * * 



0\ 



x 7 
x 3 > 



Step 3. We consider column X3 = 011 (i.e., the fourth column). Choose G^-unitary matrices U 
Uxl, ^xli ^xti ^xl with appropriate determinants to eliminate the nonzero entries in column X% 
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with row indexes X 8 , X 7 , X 6 , X 5 , X 4 . Let \J X , = Ugugugugug. Then 



Ux 3 U X2 U Xl U 
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Step 4. We consider column X4 



010 (i.e., the third column). Choose G^-unitary matrices U x 7 , 



Uy 6 , Uv 5 , f/jf * with appropriate determinants to eliminate the nonzero entries in column X4 with 



x 4 ' x 4 

row indexes X 8 , X 7 , X e , X 5 . Let Ux A 
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Step 5. We consider column X5 = 110 (i.e., the seventh column). Choose G^-unitary matrices 



U 



x 7 



X 5 



ttX 6 
U X 5 > 



U x with appropriate determinants to eliminate the nonzero entries in column X5 with 



'X 5 



row indexes Xg, X7, Xq. Let Ux 5 



uZ r >Uy 6 U% 7 . Then 

A 5 A 5 A 5 



U X5 Ux 4 Ux 3 Ux 2 U Xl U 
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Step 6. We consider column Xq = 111 (i.e., the eighth column). Choose G^-unitary U x 7 , 
with appropriate determinants to eliminate the nonzero entries in column Xq with row indexes 
X 8 ,X 7 . Let U Xe 



U% 6 UZ 7 . Then 

A6 Aj 



Ux 6 Ux 5 Ux i Ux 3 Ux 2 Ux 1 U 



/I 4 0\ 
0**00 
0**00 
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V 0001/ 



Step 7. Finally, we consider column X 7 = 101 (i.e., the sixth column). We choose a GVunitary 
Ux 7 with appropriate determinant to eliminate the nonzero entry in the (X$,X 7 ) position so that 



Ux 7 Ux 6 Ux 5 Ux 4 Ux 3 Ux 2 Ux 1 U 



(U 0\ 

10 

10 

1 

V 0001/ 



Consequently, we have U = U Xi U 



t 

x 2 



U Xr , and the construction is completed. 
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